Strain hardening in bidisperse polymer glasses: Separating the roles of chain 
orientation and interchain entanglement 
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The effects of entanglement and chain orientation on strain hardening in glassy polymers are sepa- 
rated by examining mixtures of chains with different lengths. Simulations show that the orientation 
of a molecule of a given chain length is the same in monodisperse systems and bidisperse mixtures, 
even when entangled and unentangled chains are mixed. In addition, the stress in mixtures is equal 
to the weighted average of the stresses in monodisperse systems. These results indicate that chains 
contribute independently to strain hardening, that chain orientation is determined by local inter- 
actions with the surrounding glass, and that entanglements play at most an indirect role in strain 
hardening in the range of strains typically studied. We discuss these results in the context of recent 
microscopic theories. 

PACS numbers: 61.41. +e,62.20.F-,81.40.Lm,83. 10. Rs 



I. INTRODUCTION 

The mechanical performance and failure modes of 
glassy polymers are strongly affected by strain harden- 
ing; an increase in the stress needed to deform systems as 
the strain increases. There has been great interest in un- 
derstanding the mechanisms of strain hardening so that 
it can be predicted and optimized for applications. Re- 
cent studies have challenged the traditional view that en- 
tanglements between polymer chains control strain hard- 
ening, and suggest that orientation of individual chains 
plays the dominant role. This paper uses simulations of 
polymer mixtures to separate the roles of entanglements 
and orientation. 

The most widely used 0, Q theories of strain hard- 
ening are based on rubber elasticity. They assume that 
entanglements between polymer chains act like chemical 
crosslinks in a rubber and that strain hardening arises 
from the decrease in entropy of an affinely deforming en- 
tanglement network [H- The stress is written as 



cr(A) ~ O flow + G R g(X) 



(1) 



where o~fi ow is the plastic flow stress, Gr the "strain 
hardening modulus" , A the stretch tensor describing the 
macroscopic deformation, and g a dimensionlcss func- 
tion that describes the reduction in chain entropy. This 
approach has had much success in fitting experimental 
stress strain curves, but the fits seem inconsistent with 
the underlying microscopic model pj. One difficulty is 
that Gr is two orders of magnitude larger than expected 
near the glass transition temperature T g . More troubling 
is that strain hardening grows with decreasing temper- 
ature, while any entropic force should decrease. A final 



difficulty is that the theory does not explain the plastic 
flow stress, which must be added by hand. 

A number of recent experiments and simulations have 
examined the ori gins of these difficulties [3, IE S 0, B 0, 
ES EH, El, El, Ellll, El El • One realization is that the 
two terms on the right hand side of Eq. |T]) are of simi- 
lar magnitude and must arise from similar mechanisms. 
Experiment, theory and simulations all show that ati ow 
and Gr are linearly related when pressure, T or strain 
rate is varied [H, EH, EH, EH • This connection explains 
the large magnitude of Gr and the fact that it grows 
with decreasing temperature as is typical for the plas- 
tic flow stress. A direct connection between Ufi ow and 
Gr has been established in the athermal limit (T — > 0), 
where both contributions to the stress are given by a 
single scaling factor times the rate at which interchain 
bonds break during local plastic rearrangements [H, EH • 

Simulations also show strain hardening occurs in poly- 
mers that are too short to form the entangled network 
that is assumed in rubber elasticity theories [H, EH, E3, 
E3- Moreover, the strain hardening in these unentan- 
gled systems can be mapped to that of entangled chains 
if the macroscopic deformation A in Eq. [1] is replaced 
[l2l EH by an effective stretch tensor A e / / that describes 
the orientation of individual chains: 



cr(A) = a f i ow +G° R g(X ef f), 



(2) 
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where G R is the value of Gr in the long-chain limit. The 
orientation of entangled chains is consistent with an affine 
deformation, A e // = A. For unentangled systems the 
orientation is consistent with a subaffine deformation by 

The above findings support the notion that entangle- 
ments play an indirect role and that strain hardening is 
produced by the increased rate of local plastic rearrange- 
ments that are needed to maintain chain connectivity as 
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chains orient flp . [T3L [l9j . In contrast, other recent ex- 
periments [1, |20( | suggest that entanglements play a cen- 
tral role. They find that Gr is directly proportional to 
the entanglement density p e in both pure systems and 
polymer mixtures. However, there are other explana- 
tions for this correlation. For example, an increase in 
chain stiffness increases p e and also stretches the chains 
so that more plastic rearran gem ents are required to main- 
tain chain connectivity [H O HH, Hl| • 

In this paper we examine strain hardening of polymer 
mixtures. Using polymers with different chain lengths 
but the same interactions, stiffness, etc., allows us to iso- 
late the roles of entanglements and orientation. The sim- 
ulation results show that the orientation of chains plays 
the dominant role and that the response of mixtures has 
a surprisingly simple mean-field form. The variation of 
chain orientation with macroscopic stretch is nearly inde- 
pendent of the length of surrounding chains. In addition, 
the stress in mixtures is a simple weighted average of re- 
sults for pure systems. These results suggest that strain 
hardening can be modeled as a sum of contributions from 
individual chains as they are deformed by interactions 
with the surrounding glass and produce local plastic re- 
arrangements in the glass to maintain their connectivity. 

II. POLYMER MODEL AND METHODS 

Molecular dynamics (MD) simulations are performed 
using a coarse-grained bead-spring polymer model [23[. 
All monomers have mass m and interact via the trun- 
cated and shifted Lennard- Jones potential 

U LJ {r)=A U() [{a/r) 12 -{a/rf-{a/r c ) 12 + {a/r c f] , (3) 

where r c — 1.5a is the cutoff radius and Ulj(t) = for 
r > r c . Covalent bonds between adjacent monomers on 
a chain are modeled using the FENE potential 

U(r) = -0.5fcflgln(l - ( r / R o) 2 ) , (4) 

with the canonical [23j parameter choices Rq = 1.5a and 
k = 30uo/a 2 . The equilibrium covalent bond length is 
Iq = 0.96a. We consider semiflexible chains with an an- 
gular potential Ubend — ^bend(l — cos(9)), where 9 is the 
angle between consecutive bond vectors along a chain. 
All quantities are expressed in units of the molecular 
diameter a, binding energy uq, and characteristic time 
tlj = y/ ma 2 Juq. Periodic boundary conditions are im- 
posed, with periods Lj along directions i = x, y, and 
z. 

Simulations are performed using the same protocols 
as in Refs. [ill [l2l . Il3j |. A rapid quench is applied to 
well-equilibrated melts [H| , producing glasses at the de- 
sired temperature T. The rapidity of the quench sup- 
presses aging and strain softening. The simulations pre- 
sented here are performed at T — 0.2uo/fc_B and T — 
0.275u /fcB- Both values are well below the glass transi- 
tion temperature T g ~ uo/3kB [25| . yet sufficiently high 



to observe significant thermally-activated rate-dependent 
relaxation [111 ]. Uniaxial compression is performed at 
constant strain rate e = L z /L z . The volume stays nearly 
constant during compression. In this case the function 
describing the entropy loss in Eq. [T] can be written as 
g(X) = 1/A 2 — A where A = L z jL\ is the stretch relative 
to the initial size L z q. 

The strain rates used range from 10~ 6 /tlj to 
10- 4 /Tij. While these rates are higher than experiment, 
previous studies in this range of e and T have been shown 
to capture many aspects of experiments, including log- 
arithmic rate dependence [25[, creep [2y|, and a linear 
relationship between flow and hardening modulus [l5l |. 
Our rates are also orders of magnitude slower than those 
obtainable with more computationally intensive united 
atom models (e. g. Ref. 17}). At the higher rates studied 
with these potentials, unentangled chains deform nearly 
affinely because there is insufficient time for activated 
motion of chains relative to their neighbors. This motion 
is particularly important for the chain length dependence 
of strain hardening considered here. 

Systems contain short and long chains with N s h ort 
and Ni on g monomers, respectively. The total number 
of monomers is N to t and the weight fraction of short 
chains is /. We present results for Ni ong = 350, 10 < 
Nshort < 25, and kbend = 0.75e, but equivalent results 
were obtained for a range of other values. For this kb en d, 
the entanglement length is iV e ~ 40 monomers [27|- In 
experiments, unentangled glasses with N < N e typically 
exhibit brittle fracture, but fracture is suppressed in sim- 
ulations [HI, [IH , perhaps because of the small system size 
and periodic boundary conditions. 

III. RESULTS 

Panels (a) and (b) of Figure Q] show the stress as a 
function of 5(A) for mixtures with N s h or t = 10 and 25, 
respectively. Results for / = 0.5 lie directly between the 
results for pure systems with / = and 1. Also shown is 
the weighted average of the pure system results 

&ave = (1 - f)viong + J '° 'short (5) 

for / = 0.5. Within our statistical accuracy, the stress 
in the mixtures is equal to this weighted average. We 
have verified that this "stress superposition principle" 
holds for different kbend, chain lengths, temperatures 
and /. The only discrepancies are observed at large g 
(> 3.5) when chains in entangled systems have stretched 
nearly taut and significant energy is stored in the covalent 
bonds. Experiments only reach such large stretches when 
special sample preparation procedures are used 0, H^] . 

Equation |T]) would predict a linear increase in stress 
with g. The systems in Fig.[T]were chosen to exhibit de- 
viations from this straight line behavior, but still obey 
stress superposition. Deviations from linearity in ex- 
periments are often fit by including non-Gaussian (fi- 
nite chain length) corrections to the network entropy 
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FIG. 1: Panels (a-b): Stress plotted against g(X) for uniaxial 
compression of monodisperse glasses with length Ni ong (dash- 
dotted) or N a hort (dotted) and of 50/50 mixtures of the two 
lengths (dashed). The solid curves show the average of the 
two monodisperse systems. Here kb m d = 0.75 (N e = 39). 
Panel (a): N short = 25, T = .275u /k B and |e| = W~ 5 /t L j- 
The inset shows a blowup of the data at large g [2^]. Panel 

(b) : N ahort = 10, T = 0.2u /k B and |e] = W~* /tlj. Panel 

(c) : Variation of A e // with A for pure systems with N = 350 
(solid circles), TV = 25 (solid squares), and N — 10 (solid 
triangles). Open symbols of the same shape show results for 
chains of each length in / = 0.5 mixtures. 



2\. Previous simulations are not consistent with this 
interpretation [T^j G3- They show that the upturn in 
stress for entangled systems reflects energy stored in the 
glass, while the downward curvature for unentangled sys- 
tems results from subafline deformation of the chains 
(|/n(A e //))| < |/n(A)|). Despite the different curvatures 



for the long and short chains in Fig.[T]their contributions 
to the stress in mixtures are simply additive. 

These results do not conflict with experiments that 
found Gr was proportional to p e pi [3(| • For the specific 
polymers studied, the entanglement density of mixtures 
is a weighted average of pure systems. Thus these sys- 
tems also obey stress superposition. It would be interest- 
ing to perform similar experiments on polymer mixtures 
whose entanglement density was a nonlinear function of 
weight fraction. 

The relation between p e and / is nonlinear for the 
model polymers considered here, yet plots in a previous 
paper indicated Gr was proportional to p e flTI | . While 
this is approximately true for monodisperse systems, the 
results for mixtures were misinterpreted in this earlier 
paper. A factor of (1 — /) was inadvertently excluded in 
the calculation [3l| of p e from Primitive Path Analysis 
[27l | . When p e is evaluated correctly, the hardening mod- 
ulus is no longer proportional to Pe- However, we do find 
that the values of G r from Ref. [ll| satisfy stress super- 
position. Indeed, the superposition principle explains the 
success of the previous collapse. As noted in Ref. [Til ] . 
the incorrectly determined value of p e was proportional 
to (1 — /) and so the collapse implied Gr oc (1 — /) as 
expected from superposition. The correct p e is nearly 
proportional to (1 — f) 2 and not linearly related to Gr. 

To understand why the contribution of chains to stress 
is additive, we examined the changes in molecular confor- 
mations during shear. As in Ref. the orientation was 
quantified by the effective stretch A e / / . This is defined as 
the mean stretch of chains along the compression direc- 
tion; A e // = (R z /R° z ), where R z is the rms z-component 
of the end-to-end vector of a chain and R° z the value be- 
fore strain. Panel (c) shows A e //(A) for short and long 
chains in the pure and mixed states for the same systems 
shown in panels (a-b). Remarkably, \£Ff and A^y rt do 
not themselves depend on /. In other words, the ori- 
entation of the long (short) chains is not dependent on 
whether they are in a pure system or mixed with a fairly 
high fraction of short (long) chains. This /-independence 
of A e / / holds over the same broad range of conditions as 
the stress superposition principle. 

How can we understand the apparent /-independence 
of A e //? The simplest qualitative picture is that inter- 
chain interactions in the glass tend to enforce an affinc 
deformation of the chains, but that this is opposed by 
intrachain bonds. Chains can not deform afflnely at 
small scales because this would stretch the stiff covalent 
bonds along the chains. Instead, chains are observed to 
straighten over an effective persistence length l p that in- 
creases with g [ll|, [13] • This straightening is indepen- 
dent of chain length until the associated stress becomes 
large enough to force relaxation of the chain along its 
length. For the systems considered here, the friction 
with neighboring chains is large enough that significant 
relaxation does not occur until unentangled chains are 
stretched to a substantial fraction of their fully extended 
length. Straightening chains while maintaining their con- 
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nectivity requires local rearrangements of the surround- 
ing monomers. The number of rearrangements grows 
with the persistence length and this was found to cor- 
relate directly with the strain hardening at low temper- 
atures [iM, Ea E3, H3 - 

Note that entanglements do not enter directly in this 
picture. Since entangled and unentangled chains behave 
in the same manner at small \g\, the constraints on their 
rearrangements must be determined by very local interac- 
tions with the surrounding matrix. Entanglements only 
become relevant when g is so large that l p reaches a sig- 
nificant fraction of the separation between entanglements 
[l2l [l3| . In this limit entanglements become more effec- 
tive than friction from the matrix in enforcing an afhnc 
displacement of the chain. The associated tension in the 
chains grows as l p approaches N e lo and contributes sig- 
nificant energetic terms to the strain hardening at large 

It is worth considering under which conditions this 
mean field picture should break down. For the reasons 
just described, the range of validity will be limited to 
small strains in very densely entangled systems; fortu- 
nately, by this standard, few synthetic polymers are so 
densely entangled HH. The mean field picture should 
also break down for extremely short chains where the 
high density of chain ends changes the density, friction 
and other properties significantly. Stress superposition 
will be particularly sensitive to changes in the relaxation 
rate of chains on the end-end scale, t~J , with /. These 
are known to be significant for T above T g , where T re i ax 
for short chains is increased when they are mixed with 
longer chai ns, leading to greater orientation at a given 
shear rate [33, (36| • Our results show that such coupled- 
orientation effects are small deep in the glassy state. 

IV. DISCUSSION AND CONCLUSIONS 

The simulations presented here provide further evi- 
dence that strain hardening in polymer glasses depends 
primarily on chain orientation rather than entanglement. 
They also provide new insight into the factors that con- 
trol the degree of orientation. The local coupling to the 
glassy matrix enforces an affine deformation of the chains 
(A = X e ff) a t small \g\. Hardening arises primarily from 
the plastic rearrangements needed to maintain intrachain 
connectivity as chains orient and are straightened over 
greater lengths. As the associated stress on the chains 
grows, there is an increasing tendency for them to re- 
lax along their length, leading to subaffine deformation 
(A 7^ A e //). The observation that the value of A e // for 
a given chain length is independent of / and that the 



stress due to plastic rearrangements obeys the superpo- 
sition principle, implies that chains relax nearly indepen- 
dently. The degree of relaxation depends on a competi- 
tion between the intrachain tension caused by alignment 
and friction with the glassy matrix. Entanglements only 
enter at very large \g\ where they prevent relaxation of 
chains that are much longer than N e . 

To our knowledge, no microscopic theory that predicts 
the functional form of A e / / in glasses has been pub- 
lished. Constitutive models which employ a viscoelastic 
or viscoplastic description of glassy strain hardening (e. 
g. Refs. [13, [H, typically decompose A into rubber- 
elastic and plastic parts (or use other internal state vari- 
ables), but do not explicitly account for A 7^ A e // or the 
A-dependence of nonaffine relaxation [40(. It would be 
interesting to see if a first principles theory like that of 
Ref. could be generalized to predict A e //(A). The 
qualitative ingredients that must enter such a micro- 
scopic theory are evident from the results shown here 
and in Refs. [H, [H, [lTj] • Simply put, A e // falls behind A 
when the stress in the material is insufficient to further 
affinely orient the average chain. 

Fortuitously, A e / / can now be accurately measured in 
scanning near-field optical microscopy experiments [4~0 | 
which can discriminate values for different chain lengths 
in a bidisperse system. Ref. [UJ showed A e / / falls behind 
A for entangled chains deformed slightly above T g ; similar 
studies on a variety of chain lengths below T g would be 
of great interest. These might prove that A e // is one 
of the "mesoscale" quantities [42| whose understanding 
will improve modeling of the deformation of amorphous 
materials. Other recent experiments (43| have shown that 
segmental dynamics and larger-scale polymer dynamics 
are nontrivially coupled below T g . Further developments 
along these lines would be welcome, as would extension 
of the recently improved understanding of the relation 
[44| between linear viscoelastic behavior and structural 
relaxation to the nonlinear regime. 
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